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ABSTRACT 

Previous models of dust growth in protoplanetary disks considered either uniformly laminar or 
turbulent disks. This Letter explores how dust growth occurs in a layered protoplanetary disk in 
which the magnetorotational instability generates turbulence only in the surface layers of a disk. Two 
cases are considered: a completely laminar dead zone and a dead zone in which turbulence is "stirred 
up" from the MRI acting above. It is found that dust is depleted from high altitudes in layered disks 
faster than in those cases of a uniformly laminar or turblent disks. This is a result of the accelerated 
growth of particles in the turbulent regions and their storage in the lower levels where they escape 
energetic collisions which would result in disruption. Thus the regions of a protoplanetary disk above 
a dead zone would become rapidly depleted in small dust grains, whereas the outer regions, where the 
MRI is active troughout, will maintain a small dust poplulation at all heights due to the disruptive 
collisions and vertical mixing from turbulence. This structure is similar to that which has been inferred 
for disks around TW Hydra, GM Auriga, and CoKu Tau/4, which are depleted in dust close to the 
star, but are optically thick at larger heliocentric distances. 

Subject headings: turbulence, solar system:formation, planetary systems: protoplanetary disks, 
stars:pre-main sequence 
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1. INTRODUCTION 

The initial stages of planetary system formation in- 
volve the coagulation of dust in protoplanetary disks 
into planetesimals. Models of this process have been de- 
veloped in o rder to interpret observations of protoplan- 
etary disk s ( Dullcmon d fc Dominild [20051 : ITanaka et aTl 
120051: iD'Alessio et alJl2006f > or to understand the details 
of how this process occurred in our own solar nebula 
(|Weidenschillindll980lll984 Il997f) . An important factor 
in determining how growth takes place is whether tur- 
bulence is present in the disk, as it will increase the rel- 
ative velocities between dust particles, leading to more 
frequent and energetic collisions. Previous efforts have 
focused on disks which were uniformly laminar or uni- 
formly turbulent. However, the likely source of turbu- 
lence in protoplanetary disks, the magnetorotational in- 
stability (MRI) , is expected to only operate in regions of 
the disk which are sufficiently i onized for the gas to cou - 
ple to the local magnetic field (jBalbus &: Hawlevlll991f ). 
Such regions are expected where the gas temperature 
exceeds ~1000 K allowing the gas to become collision- 
ally ionized or at the surface layers of the disk where 
absorption of X -rays and cosmic rays lead to ionization 
(|Gammidfl996l ). The region which is not ionized suffi- 
ciently, which would likely be the disk midplane in the 
planet formation region, is termed the "dead zone." 

In addition to the energetic radiation flux, the presence 
of d ust will also determ ine the ionization fraction of the 
gas. ISano et ail (|2000f ) demonstrated that grain surfaces 
allow for rapid recombination of ions and electrons, thus 
keeping the ionization fraction low. Significant ioniza- 
tion occurs only after the available surface area of grains 
is reduced, either through incorporation of grains into 
larger bodies, or by the depletion of solids by settling 
towards the midplane. 
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In this Letter, particle growth and settling in a layered 
protoplanetary disk is modeled. Two cases are consid- 
ered: that of a laminar dead zone (a—0), and that of 
a dead zone that becomes turbulent, though to a lesser 
ex tent (anz=0.1aiug;), i n agreement with the models 
of lFleming fc Stona (|2003l ). The focus here is to demon- 
strate the effect that the layered structure has on the 
vertical distribution of solids in protoplanetary disk. A 
more detailed model, which explores a wider parameter 
space will be presented in a future paper. In §2, the 
model and assumptions are described. In §3, the results 
of different simulations are presented. The implications 
for protoplanetary disk observations and structure are 
discussed in §4. 

2. MODEL DESCRIPTION 

The evolution of the number density, n, of particles 
of a given mass, m, and radius, a, at a given height, z, 
above the midplane of protoplanetary disk is given by: 
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The first term on the right hand side of the equation 
describes the vertical diffusion that takes place due to 
turbulence, where D z is the diffusivity (taken to be the 
local turbulent viscosity of the particles, acH /(l+St), 
where St is the local Stokes number of the particle), and 
p g is the local gas mass density. The second term on 
the right hand side describes the vertical settling due 
to gravity, where v z is the settling rate determined by 
the balance of gravity and the drag force due to move- 
ment through the gas. The drag force on the particles is 
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found using the laws given in iTanaka et all (2005). Fi- 
nally, the integral terms on the right hand side of the 
equation represent the changes that occur due to coag- 
ulation, where the first term describes the creation of 
particles of mass m through collisions of smaller ones, 
and the second term represents the loss of particles due 
to incorporation into larger solids. The K factors are the 
kernels, representing the product of the collisional cross 
section of the particles, 7r(ai+a2) 2 , with the relative ve- 
locities between those particles. The relative velocities 
arise due to Brownian motion, differential settling, gas 
drag (radial and azimuthal components) , and turbulence 
(when present) . In the case of tur bulence, the relative ve- 
locitie s are found as described in lDullemond fc Dominiki 
((20051) . 

Collisions may also result in disruption when velocities 
are large. A simple model is used where collisions are 
assumed to be destructive if the total kinetic energy is 
greater than the strength of the sum of the two masses, 
(rn\-\-m,2)q* . In the case of a destructive collision, all 
of the mass of the two colliding bodies is assumed to 
be dispersed into the smallest sized particles considered 
here (this is represented by the source term, S m , above). 
Here, a strength of g*=1.25 x 10 5 erg/g was assumed 
for all bodies, independent of size. This strength corre- 
sponds to bodies of equal sizes being disrupted at col- 
lisional velocities in exce ss of 10 m/s . This is likely a 
conservative estimate as IWurm et al.l (|2005[ ) found ag- 
gregation can occur at velocities as high as 25 m/s. 

The disk considered here is assumed to have sur- 
face density and temperature structures given by 
£(r)=2000(l/r AU) g/cm 2 and T(r)=300(l/r AU) 5 K 
respectively. The disk is assumed to be isothermal with 
height above the disk midplane. Inside of the snow line 
(taken to be 4 AU, where T=150 K), dust is initially 
suspended at a solid-to-gas mass ratio of 0.005 with a 
particle density of 1.5 g/cm 3 . All solids are initially dis- 
tributed as particles 0.1 /xm in diameter. In all cases, the 
disk is assumed to extend to 3 pressure scale heights, H, 
above the disk midplane. 

A simplified model for the growth of the MRI-active 
region is used here. The MRI is assumed to initially 
be active in the upper 1 g/cm 2 of a given location of 
the disk. As the available surface area of the grains de- 
creases, either through coagulation or settling, the active 
layer grows. The MRI is assumed to descend into lower 
altitudes when the total mass denity of dust smaller than 
a=l /im in a given cell is le ss than 10~ 4 time s the orig- 
inal mass density of solids (Sano et al. 2000). The ac- 
tive layer is not allowed to expand beyond 100 g/cm 2 , 
marking the maximum pen etration d e pth o f ionizing ra- 
diation, in agreement with iGammiel (|1996f ) . The MRI 
region is assumed to have a turbulence characterized by 
a=10~ 2 . In those cases specified, the dead zone develops 
a value o f a=10~ 3 when Embj to mimic the 
results of iFleming fc Stone! (|2003D . 

3. RESULTS 

Figure 1 shows the evolution of the vertical distribu- 
tion of solids at 1 AU for the cases of a uniformly laminar 
disk and a uniformly turbulent disk, which are the more 
typical cases considered when modeling dust coagulation 
in protoplanetary disks. In the case of the laminar disk, 
dust growth originally occurs due to the brownian mo- 
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Fig. 1. — Plotted are the vertical distribution of solids at 1 AU for 
a uniformly laminar disk (A) and a uniformly turbulent disk (B) at 
t=l (dotted), 1000 (dashed) 10 4 (dash-dotted), 10 5 (dash-dot-dot- 
dotted), and 10 6 (solid) years. In the case of the turbulent disk, a 
steady-state is achieved after ~1000 years, in which particles are 
created as readily as they are destroyed. 

tion that the small grains experience in the gas. As larger 
grains form, they begin to settle towards the midplane, 
with higher settling velocities at higher altitudes. These 
larger bodies "rain out" from these high altitudes, sweep- 
ing up the small grains that are present on the way to the 
midplane. As a result, the upper layers of the disk are 
first depleted in dust, with the lower altitudes becoming 
more depleted with time. Essentially all of the mass of 
the solids becomes incorporated into planetesimal-sized 
bodies which reside at the midplane. 

In the uniformly turbulent case (a=10~ 2 ), little deple- 
tion actually occurs, as turbulence inhibits the growth 
of bodies beyond a certain point. The additional 
relative velocity between particles as a result of tu- 
rublence produces a situation t hat has been referred 
to as the "meter-sized barrie r" (jWeidcnschilling 1984; 
ICuzzi fc Weidensch illing 2006) . Turbulence allows parti- 
cles near 1 meter in size to develop relative velocities with 
respect to one another that result in disruptive collisions. 
The disruptive timescale for these bodies is less than their 
growth timescale through collisions with smaller parti- 
cles. As a result, growth beyond this size is frustrated, 
and small dust is constantly resupplied to the nebula 
through these collisions and mixed vertically by turbu- 
lent diffusion. Thus, the vertical distribution of material 
reaches a steady-state situation after ~1000 years, where 
dust is somewhat depleted at high altitudes, but only by 
a factor of a few. 

Figure 2 presents the results for layered protoplan- 
etary disks, where the turbulent region begins in the 
extreme upper layers of the disk and grows with time 
as dust is depleted, with both a perfectly laminar dead 
zone and a dead zone w hich becomes turbulent when the 
IFleming fc Stonel ([2003T) criteria is satisfied. Initially, the 
dust distributions in these two cases evolve identically. In 
the MRI-active layers, particle growth is rapid compared 
to the dead zone, as turbulence increases the relative ve- 
locities between particles, leading to more frequent col- 
lisions. While growth is somewhat frustrated as larger 
bodies form and collisions become more energetic and 
disruptive, some fraction of larger bodies continuously 
settle to the deadzone, and thus are removed from the 
MRI-active layer. As there is no turbulence in the dead 
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Fig. 2. — Vertical distributions of solids at 1 AU in a layered 
protoplanetary disk for a laminar dead zone (A) and a dead zone 
which gets stirred up from the active region above (B). The lines 
represent the same times described in Figure 1. The 10 5 and 10 6 
year marks for case B are absent as all the mass of solids has been 
incorporated into large bodies at the midplane at these times. 

zone during these early times, and thus no vertical diffu- 
sion, there is no way to resupply the MRI-active region 
with solids. 

Particles in the dead zone continue to grow as described 
for the uniformly laminar disk, with much of the mass co- 
ccntrating at the midplane and being incorporated into 
larger bodies. After 14594 years, the iFleming fc Stone] 
(2003) criteria, where is met and tur- 

bulence is stirred up in the dead zone in Case B. At this 
time, much of the mass of the solids are found at the mid- 
plane in bodies 5-200 meters in radius, as seen in Figure 
3. This size range is beyond the "meter-sized" barrier. 
Thus, when the dead zone becomes turbulent, the small 
dust particles that remain can grow through mutual col- 
lisions, but only to the barrier size, at which point their 
collisions become disruptive and maintain a population 
of fine dust that gets mixed vertically. However, the tur- 
bulence also increases the rate at which dust particles 
collide with the larger bodies at the midplane, which re- 
sults in faster accretion of these particles by the larger 
bodies that have already formed. As a result, the growth 
of the large bodies is then accelerated by the presence of 
turbulence since collisions are not destructive, meaning 
that these large bodies are able to rapidly deplete the 
population of small particles that remain in the disk. 

4. DISCUSSION 

The results presented here show that in layered proto- 
planetary disks, the active regions above the dead zone 
become depleted of dust on short timescales, particu- 
larly when compared to the cases of a uniformly lami- 
nar or turbulent disk. While radial transport has been 
neglected, the timescale for dust to be delivered to the 
dead zone from above is short. Thus, it is unlikely that 
radial transport would be sufficient to replenish dust in 
the active layers to levels far above those seen here. Even 
turbulence generated by a particle rich layer at the mid- 
plane would be ineffective as the vertical scale of such 
turbulence is expected to be much less than that con- 
sidered here (fWcidcnschillingl Il997f ). This depletion of 
solids at high altitudes will have important consequences 
for modeling and interpreting observations of protoplan- 
etary disks. 
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Fig. 3. — Plotted are the mass density of particles in a layered 
disk when the MRI-active region grows massive enough to "stir up" 
turbulence in the dead zone. The lines represent the distribution 
of particles at the midplane (solid line) at one scale-height (H) 
(dashed), 2H (dash-dotted), and 3H (dash-dot-dot-dotted). 

Dust at high altitudes (3-5 scale heights) absorbs and 
scatters radiation from the central star a nd determines 
how radiation is processed by the disk (|Chiang et all 
I200H iD'Alessio etaI1l2006l ). If dust is absent from this 
region, or even depleted, it will drastica ll y affe ct this 
processing. In particular, ID'Alessio et alj (|2006l ) found 
that while dust depletions in the upper few scale heights 
would not alter the surface temperatures of the disk, mid- 
plane temperatures would decrease because less radiation 
would be processed and directed to the disk interior. This 
would be particularly important in dead zones which do 
not generate heat from viscous dissipation. 

The radial extent of the dead zone would vary from 
disk to disk, and depends on such factors as the strength 
of the magnetic field, the surface density distribution of 
the dis k, the size of the dust particles, and their compo- 
sition (jSano et al.ll2000l ). However, the general expecta- 
tion is that beyond the point where the dead zone ter- 
minates, the disk will be MRI-active at all heights, and 
therefore uniformly turbulent. Thus, this would create 
a situation where dust is depleted in the inner disk, but 
present, and at high altitudes, in the outer parts of the 
disk. Such structures have been inferred for disks around 
TW Hydra, GM Auriga, ari d CoKu Tau/4 (jCalvet et all 
120051: ID'Alessio et~aKl2005l) . These disks appear to be 
significantly depleted in dust in the inner disk, but are 
optically thick at larger radii. The transition between 
these two regions is thought to be very sharp, with the 
the inner edge, or wall, of the optically thick region being 
directly illuminated by the central star. The clearing of 
the inner parts of the disk was originally attributed to 
a planet opening a gap, but it is thought that the high 
residual mass of the disk wo uld lead to rapid in ward mi- 
gration due to disk torques (jCalvet et al.ll2005T ). 

Based on the results presented here, it is possible that 
in some disks, the inferred wall represents the termina- 
tion of the dead zone in the disk. Inside of the wall, most 
of the solids have settled to the midplane and have been 
incorporated into planctesimals, whereas outside of the 
wall, the growth of large bodies has been inhibited by 
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the presence of turbulence, as in case B of Figure 1. Not 
only would this explain why the outer regions are opti- 
cally thick, but it would also explain why dust is present 
at such high altitudes outside of the wall, but nearly ab- 
sent inside of it. The minor amount of dust inside of the 
wall in the disks around TW Hydra and GM Auriga may 
be tied to the fact that these disks are still accreting onto 
the central star. This implies that non-ne gligible mass 
trans port is taking place through the disk (| Calvet et al.l 
2005), and thus, small grains may be carried inwards 
from the optically thick outer disk by the net advective 
flows associated with disk evolution. CoKu Tau/4, on 
the other hand, is not accreting significant amounts of 
mass, and therefore would be unable to resupply the in- 
ner region with fine dust at a significant rate. For those 
disks in which no obvious inner gap appears, it is pos- 
sible that either the dead zone spans a negligible region 
of the disk, or that the dead zone is being stirred up by 
other effects. For example, gravitationally unstable disks 
have been shown to be capable of generating turbulence 
locally or to prod uce large-scale transport of fine du st on 
short timescales (|Bossll2004l: iBolev fc Durisenll2006[ ). 

These results also offer a potential way of interpret- 
ting the meteoritic record of how solids grew in our own 
solar nebula. Recent results suggest that many iron me- 
teorites are just as old as the oldest known objects in 
our sola r system, the Calci um- Aluminum- rich Inclusions 
(CAIs) (jKleine et al.ll2005f ). Previously, it was thought 
that chondritic meteorite parent bodies represented the 
oldest accreted planetesimals in our solar system. With 
chondrules, the dominant components of chondritic me- 
teorites, being typically 1-3 million years younger than 
CAIs, nebular turbulence had been invoked as a way of 
delaying the formation of large bodies until this time pe- 
riod. It has also been argued that other properties of 
chondritic meteorites, such as the narrow size range of 
chondrules in a given meteorite, are ea sier to understand 
if the y formed in a turbulent nebula (|Cuzzi "et~aLll200ll . 
120051 

The layered disk model presented here offers a possi- 
ble way to allow for the early formation of large bod- 
ies and the subsequent formation of chondritic meteorite 
parent bodies in a turbulent environment. Large bod- 
ies could form during the initial stages of disk evolution 
when the MRI operated only in the extreme upper layers 
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of the disk. The midplane region of the disk would be 
laminar and allow for the rapid formation of kilometer- 
sized bodies on relatively short timescales. These bod- 
ies could grow further through gravitational effects, and 
would likely contain enough 26 Al and 60 Fe that they 
would differentiate and develop iron cores. As the MRI- 
active region in the upper layers of the disk grew due to 
du st depletion and surpass ed the critical size as given 
by iFleming fc Stone! ([2003), the entire interior of the 
disk would become turbulent. While the case considered 
here suggested that the dust that remained in the disk 
when the dead zone became turbulent would be rapidly 
accreted by the larger bodies, this treatment neglected 
the effects of chondrule formation events, or the possi- 
ble formation of chondrite parent b odies directly from 
dust by turbulent concentration (i.e. ICuzzi et al.ll200lh . 
Additionally, if turbulence is present in the disk mid- 
plane, mass and angular momentum transport could oc- 
cur through the dead zone, allowing dust to be resupplied 
to the inner disk from larger heliocentric distances. This 
would provide a turbulent environment filled with matrix 
and chondrule precursors in which chondritic meteorite 
parent bodies could form. 

In conclusion, the upper layers of turbulently layered 
protoplanetary disks are depleted in dust more rapidly 
than either uniformly laminar or uniformly turbulent 
disks. This is because the turbulence in the upper layers 
of a disk leads to the rapid production of large bodies 
which are then stored in the laminar midplane region 
of the disk. Even if the midplane region later becomes 
turbulent because it is stirred up from above, large bod- 
ies may grow from the "seeding" that occurs during the 
initial laminar phase. While the particular results of the 
model are parameter dependent, there are aspects of pro- 
toplanetary disk observations and the meteorite record of 
our own solar system that can be understood in the con- 
text of a layered disk. Future studies should investigate 
these effects in detail. 
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